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ABSTRACT 

We report the discovery of a low-mass companion to the nearby (d = 47 pc) F7 V star HD 984. 
The companion is detected 0719 away from its host star in the L' band with the Apodizing 
Phase Plate on NaCo/VLT and was recovered by 7/-band non-coronagraphic imaging data 
taken a few days later. We confirm the companion is co-moving with the star with SINFONI 
integral field spectrograph H + K data. We present the first published data obtained with 
SINFONI in pupil-tracking mode. HD 984 has been argued to be a kinematic member of the 
30 Myr-old Columba group, and its HR diagram position is not altogether inconsistent with 
being a ZAMS star of this age. By consolidating different age indicators, including isochronal 
age, coronal X-ray emission, and stellar rotation, we independently estimate a main sequence 
age of 115 + 85 Myr (95% CL) which does not rely on this kinematic association. The mass of 
directly imaged companions are usually inferred from theoretical evolutionary tracks, which 
are highly dependent on the age of the star. Based on the age extrema, we demonstrate that 
with our photometric data alone, the companion’s mass is highly uncertain: between 33 and 
96 Mj up (0.03-0.09 M Q ) using the COND evolutionary models. We compare the companion’s 
SINFONI spectrum with field dwarf spectra to break this degeneracy. Based on the slope and 
shape of the spectrum in the //-band, we conclude that the companion is an M6.0 + 0.5 dwarf. 
The age of the system is not further constrained by the companion, as M dwarfs are poorly 
fit on low-mass evolutionary tracks. This discovery emphasizes the importance of obtaining a 
spectrum to spectral type companions around F-stars. 

Key words: stars: individual (HD984), low-mass 


1 INTRODUCTION 

Young stars are the primary targets of exoplanet imaging surveys 
because associated planets are warm and therefore bright in the in¬ 
frared. A handful of brown dwarfs and low-mass stellar compan¬ 
ions have been found in these surveys (PZ Tel B; Biller et al. 2010, 
CD-35 2722 B; Wahhaj et al. 2011, HD 1160 B,C; Nielsen et al. 
2012). Masses of directly imaged companions are estimated from 
the companion’s luminosity and theoretical evolutionary models, 
which are very sensitive to the age and distance of the host star. 

The companion mass-ratio distribution (CMRD) quantifies the 
mass ratio of a binary system (Reggiani & Meyer 2013). Based 


* Based on observations collected at the European Organization for Astro¬ 
nomical Research in the Southern Hemisphere, Chile, ESO under program 
numbers 089.C-0617(A), 089.C-0149(A), 093.C-0626(A) 


on observational data, the initial mass function (IMF) for brown 
dwarfs and very-low-mass-stars (0.08-0.2M o ) likely differs from 
stars (Thies et al. 2015). The primary formation mechanisms for 
these low-mass companions, including fragmentation and capture, 
is still under debate, making each new discovered low-mass com¬ 
panion an important test case for the theoretical formation mecha¬ 
nisms of low-mass companion. The orbital motion of the compan¬ 
ion around the primary star, measured as a small arc on the sky, 
can be used to find orbital solutions (Pearce et al. 2015). Some or¬ 
bital properties, such as eccentricity, can be constrained with only 
two measurements (Biller et al. 2010). Additionally, with multiple 
epoch astrometric measurements from direct imaging and radial ve¬ 
locity data, the dynamical mass of a companion can be measured. 
This acts as an important comparison with the inferred compan¬ 
ion masses from theoretical evolutionary models (Bonnefoy et al. 
2009; Dupuy et al. 2015; Close et al. 2007). 
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We report the detection of a companion around the F7V star 
HD 984 (HIP 1134). This star has been part of many imaging 
surveys searching for planets (Brandt et al. 2014; Rameau et al. 
2013) due to its proximity, brightness (d = 47.1 ± 1.4 pc; V = 7.3; 
ESA 1997; vanLeeuwen 2007), and proposed youth (30 Myr, 
Zuckerman et al. 2011). However, no comoving companions have 
yet been reported. Ground-based sub-mm and Spitzer infrared pho¬ 
tometry of HD 984 have not detected any evidence of a dusty debris 
disk (Mamajek et al. 2004; Carpenter et al. 2009; Ballering et al. 
2013). 

In Section 2 we describe our coronagraphic and non- 
coronagraphic observations with NaCo on the VLT. and our SIN- 
FONI/VLT integral field spectrographic data. In Section 3 we mea¬ 
sure the companion’s photometry and astrometry. In Section 4 we 
discuss previous age estimates of HD 984, estimate the age of the 
primary based on its main sequence isochronal age and other age 
indicators, and derive the mass of the companion based on evolu¬ 
tionary models and spectral analysis. We conclude in Section 5. 

2 OBSERVATIONS 
2.1 NaCo/VLT 

Observations of HD 984 were taken on UT 2012 July 18 and 20 
(089.C-0617(A), PI: Sascha Quanz) at the Very Large Telescope 
(VLT)/UT4 with NaCo (Lenzen et al. 2003; Rousset et al. 2003). 
The Apodizing Phase Plate coronagraph (APP; Kenworthy et al. 
2010; Quanz et al. 2010) was used for diffraction suppression thus 
increasing the chances of detecting a companion very close to the 
target star. Data were obtained with the L27 camera, in the L'-band 
filter (/l = 3.80/rm and AT = 0.62^m). The visible wavefront sensor 
was used with HD 984 as the natural guide star. We observed in 
pupil tracking mode (Kasper et al. 2009) for Angular Differential 
Imaging (ADI; Marois et al. 2006). We intentionally saturated the 
PSF core to increase the signal-to-noise from potential companions 
in each exposure. Unsaturated data were also obtained to calibrate 
photometry relative to the central star. 

The APP suppresses diffraction over a 180° wedge on one side 
of the target star. Excess scattered light is increased on the other 
side of the target that is not used in the data analysis. Two datasets 
were obtained with different initial position angles (P.A.) for full 
360° coverage around the target star. The field rotation was 47°.4 
in the first hemisphere and 42°.5 in the second hemisphere. 

Direct imaging observations of HD 984 were obtained on 
VLT/NaCo on UT 2012 July 20 (089.C-0149(A), PI: Julien 
Rameau). The data were taken with the L27 camera on NaCo in 
ADI pupil tracking mode. Both the saturated and unsaturated data 
were imaged in L'-band with the same exposure time, but in the 
unsaturated images a neutral density filter (ND_LONG) was used. 
The field rotation was 41°.3 in the direct imaging data. All datasets 
were obtained in cube mode. Each APP cube contains 120 frames, 
with an integration time of 0.5 s per frame. The total integration 
time in the APP was 60 min in hemisphere 1 (60 cubes) and 66 
min in hemisphere 2 (66 cubes). Unsaturated APP exposures were 
0.056 s per frame (222 frames in 12 cubes), with a total time on tar¬ 
get of 150 seconds. The direct imaging cubes contain 100 frames, 
with 0.2 s per frame, with a total integration time of 48 min for the 
saturated data and 202 seconds for the unsaturated data. 

A dither pattern on the detector was used to subtract sky back¬ 
ground and detector systematics from both datasets, as detailed 
in Ken worthy et al. (2013). Data cubes are subtracted from each 


other, centroided and averaged over. Optimized principal compo¬ 
nent analysis (PCA) is run on both of the APP hemispheres and di¬ 
rect imaging data independently, following Meshkat et al. (2014). 
Six principal components are used to model the stellar PSF, which 
results in the highest signal-to-noise detection of the companion. 
PCA processed frames are derotated and averaged to generate the 
final image with North facing up. 

2.2 SINFONI/VLT 

Data were obtained on HD 984 with the AO-fed integral field spec¬ 
trograph SINFONI (Eisenhauer et al. 2003; Bonnet et al. 2004) at 
the VLT on UT 2014 September 9 (093.C-0626, PI: G. Chauvin) 
in ~ 075 seeing. The H + K grating was used, which has a resolu¬ 
tion of ~1500. The spatial sampling was 12.5 masx25 mas (in the 
horizontal and vertical directions, respectively) resulting in a field 
of view of 078 x 078. We obtained 42 raw cubes of the target, each 
consisting of 10x4s coadded exposures. Similar to the NaCo ob¬ 
servations, data were obtained in pupil tracking (PT) mode (Hau et 
al. in prep). The companion rotated 49.26° around the center of the 
field of view during our 28 min integration. 

The instrument pipeline version 2.5.2 1 was used to correct 
the raw science frames from hot and non-linear pixels, detector 
gain, and distortion. Final cubes were reconstructed from the result¬ 
ing frames, associated wavelength map, and slitlet positions. Data 
cubes were corrected for OH lines and background emission using 
a dedicated algorithm (Davies 2007) implemented in the pipeline. 

The wavelength-dependent drift in the star position, caused by 
the atmospheric refraction, was registered (modeled by a 3rd order 
polynomial) and corrected. We took the mean of the parallactic an¬ 
gle values at the beginning and end of a given exposure, which is 
stored in the cube file headers. We processed the cubes with the 
classical-ADI (CADI) algorithm (Marois et al. 2006) to suppress 
the stellar flux independently at each wavelength (2172 indepen¬ 
dent spectral resolution elements). 

Data were corrected for telluric absorption lines using the ob¬ 
servations of an A3V standard star (HD 2811) obtained close in 
time to HD 984 and at a comparable airmass. The resulting spec¬ 
trum of HD 984 was flux-calibrated using the 2MASS Lf-band 
magnitude of the star (Cutri et al. 2003) and a spectrum of Vega. 

3 PHOTOMETRY AND ASTROMETRY OF HD 984 B 
3.1 NaCo/VLT 

The companion was clearly detected very close to the star 
(Figure 1, APP hemisphere 1 on top-left and direct imaging on 
top-right). The companion detection was confirmed using the Pyn- 
point pipeline (Arnara & Quanz 2012) and the IPAG-ADI pipeline 
(Chauvin et al. 2012). The detection was robust against changing 
the number of principal components used in the stellar PSF model. 

We used artificial negative companions to determine 
the astrometry and photometry of the companion (following 
Meshkat et al. 2015). The unsaturated stellar PSF was used to gen¬ 
erate artificial companions in the APP and direct imaging data. A 
scaling factor of 0.018 was applied to the direct imaging unsatu¬ 
rated data to account for the attenuation from the neutral density 
filter. We injected artificial negative companions into the data near 

1 http://www.eso.org/sci/software/pipelines/sinfoni/sinfoni-pipe- 
recipes.html 
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Figure 1 . Top-left: Final PCA processed image of HD 984 APP hemisphere 1 data with North facing up. Twenty principal components were used to model 
the stellar PSF in this image. Top-right: Final PCA processed image of HD 984 from direct imaging data with North facing up. Six principal components 
were used to model the stellar PSF. Bottom-left: Collapsed H + K SINFONI IFS data cubes processed with CADI, with North facing up. All three images are 
displayed in the same color scale. Bottom-right: The position of the companion is plotted as red points for both the VLT/NaCo dataset epoch (UT 2012 July 
18) and the VLT/SINFONI epoch (UT 2014 September 9). The blue point is the position of the companion if it were a background source at the time of the 
VLT/SINFONI dataset epoch (UT 2014 September 9). Error bars are included for all points. 


the expected position of the companion in steps of 0.1 pixels. The 
artificial companion contrast was varied from 5.0 to 7.0 in steps of 
0.01 mag. The chi-squared minimization over the T/D patch at the 
location of the artificial negative companion yielded the following 
results. 

Based on this analysis, the companion contrast in the APP data 
is best approximated as A L' = 6.0 ± 0.2 mag ( L' = 12.0 ± 0.2mag). 
The angular separation of the companion is 0719 ± 0.02, which 
corresponds to a projected separation of 9.0±1.0 AU for the stellar 
distance of 47.1 ± 1.4 pc from van Leeuwen (2007). The position 
angle (P.A.) of the companion is 108°.8 ± 3°.0. The error in the 
measurements is due to the range in artificial companions which 
successfully subtract the companion signal. 

The companion contrast in the direct imaging data is 
AL' = 5.9 ± 0.3 mag at 07208 ± 0.023 (9.8 ±1.1 AU). The P.A. 


of the companion is 108°.9 ± 3°.l corrected to true North ori¬ 
entation, calibrated using 6 Ori, observed at the same epoch with 
the same mode. The 6 Ori stars TCC058, 057, 054, 034, and 026 
were used to determine the true North 0.41 ±0.07° with a plate scale 
of 27.11±0.02 mas (Rameau et al. 2013). The astrometry and pho¬ 
tometry of the companion is in good agreement between the two 
datasets. Since the two datasets were only taken two days apart, we 
do not consider them separate epochs, but rather a confirmation that 
this companion is not an artifact. 


3.2 SINFONI 

The companion was detected in H + K band IFS SINFONI data 
(Figure 1 bottom-left). These data are the first published results 
demonstrating the capabilities of the SINFONI instrument in PT 
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Figure 2. Contrast curves for the SINFONI HD 984 H (black curve) and Ks 
(red curve) data, for CADI processing. The companion detection is shown 
as a point with error bars in H (black) and Ks (red), respectively. 


mode. The parallactic angle associated with each of the 43 data 
cubes was estimated by taking the mean of the parallactic angle at 
the beginning and the end of an exposure. We applied an additional 
clockwise rotation of 210.92° to the frames at the ADI reduction 
step to properly realign the field of view to the North. True North 
was estimated using GQ Lup SINFONI observations on UT 2013 
August 24, which were calibrated with NaCo observations of the 
same source on UT 2012 March 3, assuming orbital motion is neg¬ 
ligible between the epochs (following the formulae described in the 
Appendix). 

The companion is detected in all of our ADI analyses. We 
discovered that subtracting the stellar halo of the data cubes col¬ 
lapsed in wavelength followed by a realignment of the frame to the 
North also allowed detection of the companion. Unlike ADI pro¬ 
cessing, this allows for a measurement of the position of the com¬ 
panion without the biases associated with the self-subtraction of 
the companion PSF (Bonnefoy et al. 2011). We selected 10 frames 
corresponding to the first 10 cubes of the PT sequence and median- 
combined them after realigning to true North. This provided a good 
removal of residual speckles from the stellar halo. We found that 
HD 984 B lies at aPA=92.2±0.5° and a separation p = 201.6 ±0.4 
mas. The error considers uncertainties in our fitting function as 
well as true North. We also make the assumption that the instru¬ 
ment absolute orientation on sky did not vary between our UT 2014 
September 9 observations and the True North calibration (UT 2013 
August 24) with GQ Lup. The separation assumes a plate scale of 
12.5 mas/pixel reported in the instrument user manual. The plate 
scale at the time of HD 984 observations could not be measured. 

HD 984 is a relatively high proper motion star, with p a „ p s = 
102.79 ± 0.78, -66.36 ± 0.36 mas yU 1 (van Leeuwen 2007). If the 
companion were a background source, it should be due North of 
HD 984 and at a projected separation of < 071 in the SINFONI data 
(Figure 1, bottom-right). We estimate that the companion is very 
unlikely to be a stationary background source, with an estimated 


X~ probability of less than 10~ 6 based on the two epoch companion 
detections, stellar proper motion, distance estimates, and the astro- 
metric error on the SINFONI dataset (which is dependent on the 
plate scale). The new P.A. of the companion is consistent with Ke- 
plerian orbital motion. Based on the astrometry of the companion, 
we confirm it is bound to the star and not a background object. 

The companion flux was integrated in a 6 pixel wide circu¬ 
lar aperture in the processed CADI cubes to generate the spectrum. 
We applied the same procedure to HD 984. We corrected for flux 
losses in the companion’s spectrum (caused by the image process¬ 
ing algorithms) by adding artificial sources. We used the unsatu¬ 
rated primary star itself to scale and inject artificial sources in the 
data. Artificial sources were added in each wavelength at the same 
projected separation as the companion, but with a P.A. difference 
of -90, +90, and +180°. The resulting spectrum was corrected for 
the telluric lines using the primary spectrum. We obtained the fi¬ 
nal flux-calibrated spectrum of the companion by multiplying the 
flux ratio between the system components and the flux-calibrated 
spectrum of the star. Based on the flux ratio between HD 984 and 
B in H and AA-band, we estimate H 2 mass = 12.58 ± 0.05 mag and 
Ks 2 mass - 12.19 ± 0.04 mag for the companion. The spectrum of 
the companion is analyzed in Section 4.2. 

Figure 2 shows contrast curves comparing the sensitivity to 
point sources in the data processed with the CADI algorithm 
in H and AA-band. The contrast curve was generated following 
Chauvin et al. (2015). We injected fake planets every 10 pixels radi¬ 
ally at P.A.s of 0, 120, and 240° between 125 and 375 mas, in order 
to correct for flux losses. The fake planets are created by scaling 
the flux of the primary star. We created a pixel-to-pixel noise map 
by sliding a box of 5 x 5 pixels from the star to the limit of the SIN¬ 
FONI field of view. The 5cr detection limit is found by dividing the 
pixel-to-pixel noise map by the flux loss, taking into account the 
relative calibration between the fake planet and the primary star. 

4 ANALYSIS 

4.1 Primary Star HD 984 

HD 984 is a V = 7.32 mag (ESA 1997) F7V star (Houk & Swift 
1999) at a distance of 47.1 ± 1.4 pc (van Leeuwen 2007, m = 
21.21 ±0.64 mas). The star has color B-V = 0.522 ±0.010 from the 
Tycho catalog (ESA 1997), and for its calculated absolute magni¬ 
tude (My =3.95 ± 0.07) the star lies squarely on the field star main 
sequence of Wright et al. (2004), which predicts My(MS) = 3.95 
for this B-V color. It is not listed in latest version 2 of the Washing¬ 
ton Double Star catalog (Mason et al. 2001 ). The interstellar red¬ 
dening and extinction towards HD 984 is likely to be negligible as 
it lies within the Local Bubble (e.g. Reis et al. 2011). Examining 
the reddening estimates for 19 stars in the catalog of Reis et al. 
(2011) that have Galactic longitude within ±20° of HD 984 and 
latitude b < -48° (their survey only extends to b = -60°), we find 
no stars within 60 pc that have a E(b-y) reddening that exceeds 1 <x 
from zero. Indeed, a fit of distance versus E(b-y) reddening for the 
19 stars in this region of sky with d < 150 pc is consistent with 
negligible reddening (E(b-y) = d pc x 0.01 mmag). Given this trend, 
for a star at d — 47 pc, one would expect reddening of E(b-y) = 2 
mmag and extinction A v = 0.01, i.e. utterly negligible compared to 
our photometric uncertainties. Based on this, we assume negligible 
reddening or extinction for HD 984 and its companion. 

2 8 June 2015 version at http://vizier.cfa.harvard.edu/viz-bin/Cat7B/wds. 
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Numerous estimates of the effective temperature have 
been reported (e.g. Valenti & Fischer 2005; Masana et al. 2006; 
Schroder et al. 2009; Casagrande et al. 2011). We adopt the recent 
value from Casagrande et al. (2011) (6315 ± 89 K), which is also 
close to the median value among recent published estimates. From 
fitting the optical-infrared photometry to the dwarf color sequences 
from Pecaut & Mamajek (2013), we derive an apparent bolomet- 
ric magnitude of FID 984 to be m bo i = 7.257 ± 0.01 and bolometric 
flux 3 of 31.94 ± 0.29 pWm~ 2 . Employing the van Leeuwen (2007) 
parallax, this bolometric flux translates to a luminosity of log(L/L 0 ) 
= 0.346 ±0.027 dex. Adopting the Casagrande et al. (2011) T e ff, 
this is consistent with a stellar radius of 1.247 ± 0.053 Rq. 

4.1.1 Previous Age Estimates 

Zuckerman et al. (2011), Malo et al. (2013), and Brandt et al. 
(2014) consider FID 984 to be a member of the Columba group 
with an age of 30 Myr (Torres et al. 2008). Malo et al. (2013) 
has demonstrated that FID 984 appears to be comoving with the 
Columba group. Flowever, there remains the possibility that FID 
984 could be a kinematic interloper or that the Columba group 
is not sufficiently characterized to reliably assign ages based on 
kinematic membership. Isochronal ages of <0.48 Gyr (68%CL; 
Takeda et al. 2007), 1.2Gyr (Valenti & Fischer 2005), and 
3.1+} ° Gyr (Holmberg et al. 2009) have been estimated. 

4.1.2 Isochronal Age 

The HR diagram position of HD 984 is plotted in Figure 3 along 
with isochrones from Bressan et al. (2012) (for Z = 0.017). HD 984 
is a main sequence star, and any derived isochronal age will have 
large uncertainties. Simulating the HR diagram position and inter¬ 
polating their ages and masses, the isochronal ages are consistent 
with 2.01};} (95%CL) Gyr and 1.20 ±0.06 (95%CL) M 0 (dashed 
line). However, as Figure 3 shows, a 2cr deviation (dashed error 
bar line) in both T c ff and log(L/L e ) are consistent with the 30 Myr 
isochrone. It takes a 1.2 M 0 star 27 Myr to reach the main sequence, 
and ages of <30 Myr can be ruled out based on the isochronal age 
constraints for the secondary HD 984 B (see Section 4.2). Based on 
the HR diagram position of HD 984 compared to stellar evolution¬ 
ary tracks and isochrones, we deduce that the isochronal age of HD 
984 is likely between ~30 Myr (ZAMS) and ~4 Gyr. 

4.1.3 Other Age Indicators 

We examine other age diagnostics for HD 984 to help constrain 
its age. HD 984 appears to be fairly chromospheric ally and coro- 
nally active, which hint at fast rotation and youth. Wright et al. 
(2004) estimated an age for HD 984 of 0.49 Gyr based on chro¬ 
mospheric activity. The star has mass ~1.2 M 0 , B-V = 0.52, and 
T eff ^ 6300 K, which corresponds almost exactly with the “Kraft 
break” (Angus et al. 2015; Kraft 1967). This region divides F-type 
stars with thin convective envelopes, mostly radiative interiors with 
small convective cores, and which slow down less efficiently dur¬ 
ing their main sequence phases, from solar-type stars with thicker 
convective envelopes, radiative cores, and which slow down ap¬ 
preciably during their main sequence phase. The gyrochronol- 
ogy relations are not well-constrained near the Kraft break (e.g. 

3 Our F to! value is within 3.5% of that of Casagrande et al. (2011): 30.824 
pW m 2 . 



Figure 3. HR diagram position of HD 984 and the Sun, with isochrones and 
evolutionary U'acks from Bressan et al. (2012) overlaid (z = 0.017), with fir 
solid line and 2<x dotted line error bars. The solid line is a 30 Myr isochrone 
corresponding to the purported age of the Columba group. Dotted lines are 
isochrones with ages of 1, 2, 3. 4. 5 Gyr. The dashed line is a 1.20 ±0.06 
(95%CL) M 0 tack. 


Mamajek & Hillenbrand 2008) - and may not be applicable. We 
can, however, quantify just how active HD 984 is compared to 
nearby stars of similar spectral type. 

Coronal X-ray and Chromospheric Ca H &- K Emission: 
HD 984 has an X-ray counterpart (1RXS J001410.1-071200) in 
the ROSAT All-Sky Survey Bright Source Catalog (Voges et al. 
1999). The X-ray source is 4” away from the optical position, 
but with position error of 8”, the match with HD 984 is highly 
probable. Based on the ROSAT X-ray count rate (0.265 ± 0.029 
cts~') and hardness ratio HR1 (HR1 = -0.20 ±0.10), and using 
the formulae in Fleming et al. (1995) and Mamajek & Hillenbrand 
(2008), we estimate the following properties: 0.1-2.4 keV X-ray 
flux fx = 1.92 x 10- 12 erg s _1 cm -2 , log(L Y /erg/s) = 29.71 dex, and 
log(L Y /L to; ) = -4.22 dex. It is unlikely that the X-ray emission is 
from the low-mass companion alone as the star would have to be 
emitting ~40x more powerful than that predicted for saturated X- 
ray emission (log(L Y /L AD/ ), so the primary is almost certainly com¬ 
pletely dominating the X-ray flux. Among nearby F-type dwarfs 
in the ROSAT study of Hiinsch et al. (1999), the mean relation 
of coronal X-ray activity values goes as logfL^/L;,,,;) ^ -6.60 + 
2.67(B-V), withrms scatter ±0.68 dex. HD 984 is 0.98 dex (~1.4cr) 
brighter in log(L Y /Lt 0 /) than this locus and suggestive that the star 
is amongst the most active ~4% of F-type dwarfs. 

Based just on the coronal X-ray emission {\og{Lx[C bo i = -4.22 
dex), one would predict strong chromospheric activity logfRJ^) 
^ -4.34 using equation Al of Mamajek & Hillenbrand (2008). 
Indeed, this predicted log(R^ A ,) is close to that reported from 
multiple surveys - logfRj^) ranges from -4.31 (Schroder et al. 
2009) to -4.46 (Pace 2013), with a mean value of approximately 
-4.40 (Wright et al. 2004; White et al. 2007; Schroder et al. 2009; 
Isaacson & Fischer 2010; Pace 2013). In a catalog of unique FGK 
stars with logfR^) values from the surveys of Henry et al. (1996), 
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Wright et al. (2004), and Gray et al. (2006), one finds 146 sin¬ 
gle main sequence Hipparcos stars with B-V within ±0.02 mag 
of HD 984’s color (0.52). Among these dwarf stars of similar 
color to HD 984, its loglRj^) value (~ -4.40) places it among 
the most chromospherically active 5%, similar to that inferred 
from the coronal activity. The star’s combination of color (B-V = 
0.52) and chromospheric activity (logfR^) =; -4.40) for HD 984 
can also be compared to stars to various samples in Fig. 4 of 
Mamajek & Hillenbrand (2008): Sco-Cen (~0.01-0.02 Gyr), the 
Pleiades (—0.13 Gyr), the Hyades (~0.6 Gyr), the M67 cluster (~4 
Gyr), and field stars. The star’s chromospheric activity is less active 
than typical Sco-Cen members, most consistent with the Pleiades 
locus, would be near the upper limit of that seen among Hyades 
members, and is way more active than M67 members and the me¬ 
dian values calculated for late F-type stars. 

Rotation: Given the star’s projected rotational velocity (usin; 
= 42.13 ± 1.65 kins -1 ; White et al. 2007) and the calculated ra¬ 
dius (1.247 ± 0.053 R g ), this implies that the star’s rotation period 
must be P ro , <1.5 day. Note that among 487 stars classified as F7 
dwarfs in the Glgbocki & Gnacinski (2005) compendium of pro¬ 
jected rotational velocities (usin;), HD 984 is among the top ~1% 
fastest rotating stars. The F7 dwarfs in the Glgbocki & Gnacinski 
(2005) catalog have median usin; ^ 7.2 kms -1 with 95% range 
of 2.7-35.2 kms -1 . The mean usin; value for a F7 star of simi¬ 
lar color and absolute magnitude to HD 984 in the ~0.6 Gyr-old 
Hyades cluster is approximately 12kms -1 , and none within ±0.05 
mag in B-V of HD 984 have usin; > 20 km s -1 (Soderblom et al. 
1993). In the Pleiades, however, there are stars of similar color 
to HD984 with usin; of ~40-60 kms -1 (Soderblom et al. 1993). 
The star’s combination of log(L x /Lu D ;) and usin; for its V-I color 
(0.59; ESA 1997) appears typical among F-type stars in the ~50 
Myr-old IC 2391 and IC 2602 clusters, and would place it among 
the fastest rotators in the ~130 Myr-old Pleiades (see Fig. 6 of 
Stauffer et al. 1997), where the quoted ages are Li-depletion ages 
from Barrado y Navascues et al. (2004) and Dobbie et al. (2010). 

The rotation, coronal activity, and chromospheric activity in¬ 
dicators for HD 984 seem consistent with a very young, active main 
sequence star. The rotation and X-ray emission are consistent with 
the star being amongst the youngest ~l-5% of F7 dwarfs. The 
main sequence lifetimes of ~1.2 M 0 stars are ~5 Gyr (Bressan et al. 
2012), suggesting that if HD 984 is amongst the youngest ~l-4% 
of F7 dwarfs, that its age is very likely to be <200 Myr. This would 
be consistent with its chromospheric activity being similar to that 
of Pleiades members, its rotation being similar to that of ~50 Myr 
stars in IC 2391 and IC 2602 (although near the upper limit of pro¬ 
jected rotation velocities for Pleiades members), and its rotation 
and activity indicators all being very inconsistent with Hyades- 
age stars. Based on rotation and activity indicators, we adopt a 
strong upper limit on the age of HD 984 of <200 Myr. The fact 
that the star is on the main sequence places an approximate lower 
limit of 30 Myr. We conclude that an approximate age for the pri¬ 
mary star HD 984 that is independent of any kinematic membership 
or consideration of the properties of the secondary is 30-200 Myr 
(95%CL). 

4.2 Companion Characteristics 

Based on the NaCo L' photometry alone, the companion mass can 
only be estimated from theoretical models which are highly depen¬ 
dent on the age of the star and lead to large uncertainties in the de¬ 
termined mass. Using the COND evolutionary tracks (Baraffe et al. 
2003), if the system is 30 Myr, the companion may be a 33 ± 6Af Jup 
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Figure 4. Comparison of HD 984 B spectrum (black) to those of field 
dwarfs, young companions, isolated objects, and to the best-fitting BT- 
COND spectrum. It enables us to conclude that HD 984 B is an M6.0 ± 0.5 
dwarf. 


young brown dwarf. If the system is 200 Myr, the companion may 
be a 0.09 ± 0.01 M e low mass star, likely an M5-M6 dwarf. The 
errors are based on the uncertainty in the photometry and the dis¬ 
tance of the star and do not include systematic uncertainties in the 
models. The range of companion masses based on the DUSTY evo¬ 
lutionary model (Chabrier et al. 2000) is within the errors of the 
values derived from the COND model: 30 Myr and 34 ± 6Mj up or 
200 Myr and 0.10± 0.01M o . Thus, the companion mass estimation 
based on photometry is not significantly impacted by the evolution¬ 
ary model chosen. 

We compare our CADI analysis of the SINFONI spectrum 
(black spectrum, see Figure 4) with field dwarfs from the IRTF li¬ 
brary (Rayner et al. 2009). We perform a least squares fit of the 
spectra to determine the best fit. In addition to the overall spec¬ 
tral slope from 1.5 to 2.45 fj m, we aim to fit the following spec¬ 
tral features: the Ki band at 1.516 yum, the Cai triplet near 1.98 
gm, the Nai doublet near 2.207 gm, the CO overtones longward 
of 2.3 yum, and the overall shape of the AA-band which is sensitive 
to the collision-induced absorption of molecular hydrogen, thus to 
the atmospheric pressure and surface gravity. The companion has 
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Figure 5. HR diagram position of HD 984 B in bolometric lumi¬ 
nosity (right) based on our spectral type estimate and adopting the 
Stephens et al. (2009) T c ff scale (circle), BT-COND and GAIA syn¬ 
thetic spectra (Brott & Hauschildt (2005); Allard et al. (2013): square). The 
BHAC15 (Baraffe et al. 2015) evolutionary tracks are plotted for 30, 50, 
100, 200 Myr and 1 Gyr. The 30 and 200 Myr isochrones are in blue to 
indicate the consistency with the lower and upper age limits determined in 
this work. 


Table 1. System properties 


Property 

HD 984 


HD 984 B 

Distance (pc) a 


47.1 ±1.4 


Age (Myr) b 


30-200 


Av c 


0 


Teff 

6315 ±89 


9777+127 d 

' -130 




2900 ± 200 e 

Spectral type 

F7V 


M6.0 ± 0.5 

log (L/Lo) 

0.36 ±0.03 


-2.815±0.024 

Mass 

1.20 ±0.06 Mq 


74 + |[ Mj up d 
90 + | M Jup c 
84j® Mj up f 

Separation (mas) 


201.6 ±0.4 

PA.(°) 

H 

6.170 ±0.023 

92.2 ± 0.5 

12.58 ±0.05 

Ks 

6.073 ±0.038 


12.19 ±0.04 

V 

6.0 ±0.1 


12.0 ±0.2 


a Hipparcos catalog (van Leeuwen 2007). 
b This work. 

c Schlegel et al. (1998) and Reis et al. (201 1). 

<l Based on the spectral type conversion scale of Stephens et al. (2009). 
e Based on BT-COND and GAIA synthetic spectra. 
f Converted from the derived log(L/L 0 ) for HD 984 B. 


features later than M5 and midway between those of M5.5 and 
M6.5 field dwarfs. The SINFONI spectra of the young (<150 Myr) 
mid-M companions AB Dor C and TWA 22 B (Close et al. 2007; 
Bonnefoy et al. 2009) reproduce the pseudo-continuum in the H- 
band, but have bluer slopes than our companion, suggesting a later 
spectral type for HD 984 B. Spectra of other young mid-M dwarfs 
have a more triangular //-band shape than our object. The com¬ 
panion fits well with the M6 object SERC 296A (Thackrah et al. 
1997; Allers & Liu 2013), which does not appear to be a good can¬ 
didate member to any nearby associations based on its kinematics 
(Gagne et al. 2014), but has an age below ~ 200 Myr due to lithium 
absorption. We therefore conclude that the HD 984 B is more likely 
a M6.0 ± 0.5 object, which is younger than the typical field dwarf 
age (>>1 Gyr). The subtype accuracy is due to the companion fea¬ 
tures being intermediate between a M5.5 and an M6.5. 

The spectral type of M6.0±0.5 corresponds to T cff = 2777*}^ 
K using the conversion scale of Stephens et al. (2009). The T c ff ac¬ 
curacy is the quadratic combination of the spectral type estimate 
and systematic uncertainty. Conversely, the companion’s pseudo¬ 
continuum shape and main absorption features are best reproduced 
by BT-COND and GAIA synthetic spectra (Brott & Hauschildt 
2005; Allard et al. 2013) with T cff = 2900 ± 200 K and log g=5.0- 
5.5 dex. The accuracy on T e ff is limited by the intrinsic differences 
and uncertainties between the atmospheric model grids (//-band 
is known to be badly reproduced by models, see Bonnefoy et al. 
2013). Rajpurohit et al. (2013) find a T e ff range of 2700 to 3000 
K for M5.5 to M6.5 dwarfs by fitting BT-Settl (Allard et al. 2013) 
spectra to optical spectra of M-dwarfs. This is consistent with our 
spectral type and derived temperature. 

We calculate the bolometric luminosity based on our de¬ 
rived spectral type and H, Ks and L'-band companion magni¬ 
tude measurements. We adopt bolometric corrections of BC H (B) = 
2.588 ± 0.032 mag, BCk S (B) = 2.940 ± -0.015 mag, and BC l ^b) = 
3.260 ± 0.022 mag by interpolating between an M6 and M7 
dwarf, using bolometric correction and color information for M 
dwarfs from Pecaut & Mamajek (2013), Schmidt et al. (2014), and 
Dupuy & Liu (2012). By taking the weighted mean of our bolomet¬ 
ric luminosity calculations from H, Ks and L'-band measurements, 
we find log(L/L Q )=-2.815±0.024 dex. These results are summa¬ 
rized in Table 1. 

Figure 5 compares the derived T c ff of the companion against 
the BHAC15 (Baraffe et al. 2015) evolutionary tracks of different 
ages. Comparing our derived bolometric luminosity with evolu¬ 
tionary tracks suggests that the companion is consistent with our 
derived stellar age of 115 ±85 Myr. HD 984 B demonstrates the 
challenges in fitting M-dwarfs on low mass evolutionary tracks. 
Thus, we can rule out ages less than 30 Myr based on the com¬ 
panion HR diagram position and we have demonstrated that this 
companion is consistent with an M6.0 ± 0.5 object. 

Recently, Hinkley et al. (2015) announced the detection of 
seven companions with short projected separations around inter¬ 
mediate mass stars in the Scorpius-Centaurus association. HD 984 
B has a companion-to-star mass ratio of 0.067^'^ based on the 
mass derived from the luminosity. Given the small separation of 
the system (9.5 AU), HD 984 B shares similar characteristics as the 
aforementioned companions discovered by Hinkley et al. (2015). 
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5 CONCLUSION 

We report the discovery of a low-mass companion to the F7V star 
HD 984. This companion was detected in L'-band with the Apodiz- 
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ing Phase Plate coronagraph and non-coronagraphic photometry 
with the NaCo instrument on the VLT. HD 984 has been reported 
to be part of the 30 Myr old Columba association. However, based 
on our independent analysis of several age indicators, we estimate 
a main sequence age of 115 ± 85 Myr (95% CL). Due to the age un¬ 
certainty, the companion mass may range from a low mass brown 
dwarf (—33 M Jup ) to an M dwarf (~0.09 M 0 ), using the COND evo¬ 
lutionary models. 

We analyze the slope and shape of SINFONI H + K IFU data 
of the companion compared with field M dwarfs. We conclude the 
companion is an M6.0 ± 0.5 dwarf. Using our derived spectral type, 
we aim to determine the age of the system by placing the com¬ 
panion L'-band absolute magnitude and bolometric luminosity on 
COND evolutionary tracks. While the L'-band HR diagram posi¬ 
tion allows us to rule out an age less than 50 Myr, the companion’s 
bolometric luminosity position is consistent with an age of >30 
Myr. Thus, we cannot set age constraints on the companion, due to 
this discrepancy between the companion’s position on evolutionary 
tracks. Future observations in ./-band could help to look for signa¬ 
tures of low surface gravity in HD 984 B. 

Given its small projected separation of just ~9 AU, HD 984 
B will show significant orbital motion over the next few years, al¬ 
lowing the potential for dynamical mass determination. In order to 
determine the individual mass components of the HD 984 system, 
RV data must be obtained. Using the current projected separation as 
an approximation for the semi-major axis and assuming a circular 
orbit, we would expect to detect a maximum semi-amplitude from 
RV of 1.01 km/s for an M-dwarf companion and 0.27 km/s for a 
brown dwarf. This calculation assumes the orbit is edge-on (z'=90°), 
however based on our two epoch astrometric measurements, its or¬ 
bit is unlikely to be edge-on. Even if the companion is closer to 
face-on (z=5°), the semi-amplitude of 0.087 km/s for an M-dwarf 
and 0.024 km/s for a brown dwarf is within the detection limits of 
RV. This system is also a promising target for GAIA, which could 
reveal the reflex motion of the star and assess the dynamical mass 
of the companion. These measurements are amplified if the system 
is face-on, complementing the RV method. Thus, the dynamical 
mass of this companion can be achieved with RV and astrometric 
measurements, providing a crucial comparison with the theoretical 
evolutionary models for mass determination. 

We have demonstrated that, given the difficulty in deriving 
a reliable age for HD 984 that is independent of its purported 
group membership, the derived color and spectral parameters of 
the companion are necessary to determine the companion mass. 
These results suggest that caution should be used when estimating 
the masses of companions based on photometric data and stellar 
age based on kinematic group membership alone. It reinforces the 
importance of future near-infrared high contrast integral field spec¬ 
trographs to the characterization of low-mass stellar and substellar 
companions. 
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APPENDIX A: SINFONI ORIENTATION 

We calibrated the absolute orientation of the field of view of SIN¬ 
FONI using observations of GQ Lup B from UT 2013 August 24 
(technical program ID 60.A-9800). These observations were ob¬ 
tained using the pre-optics offering a 50x100 mas sampling in the 
H + K band. During the 54 min pupil tracking (PT) sequence, 
100x7s integrations were recorded. The field rotated by 12.84°. 
We reduced these data following the same procedure as HD 984 B. 
Once the cubes were corrected for atmospheric refraction, we re¬ 
moved the halo from the primary star centered in the field of view 
with a radial profile. 

We found that a clockwise rotation by the ADA.POS ANG + C 
values of each data cubes re-aligned the final GQ Lup data cubes 
with the North. ADA.POS ANG is a variable found in the image 
header, corresponding to the position angle of the rotator at the 
Cassegrain focus at the time of the observations. C is an additional 
offset related to the calibration of the instrument rotator true North 
position. 

We define ROT.PT.OFF + C = 180° + ADA.PUPILPOS + C, 
the angular offset needed to realign the frames to the North when 
the parallactic angle is 0. ADA.PU PILPOS is a keyword stored 
into the file header, which depends on the telescope pointing po¬ 
sition and time of observation. It is redefined at the beginning of 
any PT sequence, but it remains constant during a PT observing 
sequence. 

We verified that this relation remains valid for other datasets 
obtained in PT mode on GQ Lup during the same night, on AB Dor 
C (UT 2013 October 17, 60.A-9800), and on HD 984 B (UT 2014 
October 10, 2014 December 3,5,8 094.C-0719), and of the astro- 
metric binary HD 179058 AB (UT 2014 April 26, 60.A-9800). We 
note that we could not use the observations of the HD179058 AB 
to properly calibrate the instrument plate scale and absolute orien¬ 
tation since the binary had likely moved on its orbit since its latest 
independent astrometric measurement (Tokovinin et al. 2010). 

For the case of GQ Lup B, the ADA.PU PILPOS keyword was 
fixed to -0.07297. Therefore, we adopted ROT.PT.OFF= 179.927. 
We estimated a value of C=0.0 ± 0.5° comparing the resulting po¬ 
sition angle in the derotated cubes to the position angle of the sys¬ 
tem measured from VLT/NaCo data obtained on 2012 March 3 and 
reported in Ginski et al. (2014). We note that this value of C as¬ 
sumes that the companion did not have significant orbital motion 
in the course of one year. This is reasonable given the available 
VLT/NaCo astrometry of the system recorded since 2008 (see Ta¬ 
ble 2 of Ginski et al. 2014). We also deduce from the SINFONI 
data of GQ Lup B that the mean square plate scale is 49.30 ± 0.14 
mas/spaxel when the 50x100 mas and H + K band mode of the 
instrument are chosen. 
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